The Goldstone mode of planar optical parametric oscillators 
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We propose an experimental setup to probe the low-lying excitation modes of a parametrically 
oscillating planar cavity, in particular the soft Goldstone mode which appears as a consequence 
of the spontaneously broken (7(1) symmetry of signal/idler phase rotations. Differently from the 
case of thermodynamical equilibrium, the Goldstone mode of such a driven-dissipative system is an 
overdamped mode, whose linewidth goes to zero in the long-wavelength limit. When the phase of 
the signal/idler emission is pinned by an additional laser beam in the vicinity of the signal emission, 
the (7(1) symmetry is explicitely broken and a gap opens in the Goldstone mode dispersion. This 
results in a dramatic broadening of the response to the probe. Quantitative predictions are given 
for the case of semiconductor planar cavities in the strong exciton-photon coupling regime. 

PACS numbers: 89.75.Kd, 11.30.Qc, 42.65.Yj, 71.36.+C, 
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I. INTRODUCTION 



A central concept of the modern theory of phase tran- 
sitions in systems at thermodynamical equilibrium is the 
so-called Goldstone mode, which appears as a conse- 
quence of the spontaneous breaking of a continuous sym- 
metry Unless the system has long-range interactions 
, this mode has the general property of being a soft- 
mode, whose frequency dispersion wc(fc) tends to zero in 
the long-wavelength limit k — ► 0. The Goldstone mode 
corresponds in fact to a spatially slow twist of the or- 
der parameter, whose energy cost tends to zero in the 
long- wavelength limit. Among the most celebrated exam- 
ples of Goldstone modes in condensed matter physics we 
can mention the zero-sound mode of supcrfluid Helium 4 
or dilute Bose-Einstein condensates and the magnon ex- 
citations in ferromagnets. Zero-sound is related to the 
spontaneous breaking of the U(l) gauge symmetry of 
the quantum Bose field below the Bose-Einstein conden- 
sation temperature 0, 0, Q, while the magnon branch 
is related to the spontaneous breaking of the rotational 
symmetry of the magnetic moment orientation below the 
Curie temperature |5|. 

The concept of a Goldstone mode plays an important 
role also in the physics of nonlinear dynamical systems far 
from thermodynamical equilibrium [6j, whose stationary 
state is not determined by a thermodynamical equilib- 
rium condition, but is rather the result of a dynamical 
equilibrium between an external driving force and the 
dissipation. In many examples of such driven-dissipative 
systems a homogeneous state develops a non-trivial spa- 
tiotemporal pattern when the driving force exceeds a 
critical value. The most famous example of this pat- 
tern formation behaviour is perhaps the regular periodic 
arrangement of Benard cells that appears in heat con- 
vection through a viscous fluid when a sufficiently large 
temperature gradient (driving force) is applied in the ver- 
tical direction so to exceed the braking effect of viscosity 
(dissipation) Q • The continuous translational symmetry 
of the initially homogeneous system is reduced to the dis- 



crete symmetry of the periodic pattern of Benard cells. 
As a consequence of the spontaneously broken symmetry, 
a neutral mode of vanishing frequency and damping rate 
appears in the linear stability analysis around the sta- 
tionary state of the system, a mode which corresponds 
to the rigid translation of the whole roll pattern in space. 
This neutral mode is the non-equilibrium counterpart of 
the Goldstone mode of equilibrium statistical mechanics. 

Another well celebrated example of pattern formation 
takes place in optical parametric oscillation (OPO) in 
planar cavities ||. Differently from the more usual case 
of spherical mirror cavities with a discrete and well- 
spaced set of modes, planar cavities dispose of a con- 
tinuum of modes in which the parametric emission can 
take place, so that the light field has a rich spatial dy- 
namics 9J. In the OPO state, a periodic spatial pattern 
is indeed created in the cavity plane by interference of the 
pump, signal and idler fields. In addition to their interest 
from the point of view of fundamental physics, technolog- 
ical applications of planar OPOs have also been actively 
investigated in view of using them as flexible light sources 
in new frequency ranges, or even as sources of entangled 
photons for quantum cryptography |lfj| . 

In recent years, a great deal of theoretical and ex- 
perimental activity has concerned planar semiconduc- 
tor microcavities with quantum well excitonic transitions 
strongly coupled to the cavity mode 0, 0] . At linear 
regime, the elementary excitations of these systems are 
cavity-polaritons, i.e. a superposition of a cavity photon 
and a quantum well exciton. Cavity-polaritons have the 
interesting property of combining the extremely strong 
Kerr nonlinearity due to the excitonic component to a 
peculiar dispersion relation which allows for easy and ro- 
bust phase-matching of the optical parametric process. 
In this way, large values ofgain have been observed in 
parametric amplification [13j, as well as parametric os- 
cillation at low pump intensity values 0, 0|. Exten- 
sive studies of the parametric oscillation process in semi- 
conductor planar microcavities have verified the coherent 
nature of the parametric emission above threshold by ob- 
serving a frequency narrowing of the emission |l4j |. the 
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absence of dispersion , as well as the long-range spa- 
tial coherence 16] . 

As it happened for atomic Bose-Einstein conden- 
sates , a good deal of information on the state of the 
polariton system can be obtained by looking at the ele- 
mentary excitations of the system around its dynamical 
equilibrium state. The elementary excitations around the 
pump-only state below the threshold for parametric os- 
cillation have been theoretically investigated in , and 
predictions have been put forward for polariton superflu- 
idity effects. In the present paper, we shall investigate 
the elementary excitation spectrum around the paramet- 
ric oscillation steady state above threshold. Since a con- 
tinuous U(l) symmetry related to the signal/idler phases 
is spontaneously broken, a Goldstone mode has to appear 
in the excitation spectrum of the planar system, whose 
frequency and damping rates go to zero in the k — > 
long- wavelength limit. Although its existence has been 
recently mentioned 0] and its role in the destruction of 
long-range coherence in ID systems discussed 0, |2jj, 
no experiments have been so far performed nor proposed 
which are able to directly observe the Goldstone mode 
of optical parametric oscillators. Pioneering work on the 
elementary excitation spectrum around a parametrically 
oscillating state was reported in Ref.[2l| both in the pres- 
ence of a laser beam driving the signal mode, and in a 
pump-only configuration. The emission was experimen- 
tally observed and its spectrum compared to a calcula- 
tion of the elementary excitation dispersion. However, 
no specific attention was paid in the experiment to the 
region where the Goldstone mode is expected to appear, 
nor any mention made to it in the theoretical analysis. 
The same with other recent theoretical papers that ad- 
dress the excitation spectrum around a parametri cally 
oscillating state of a semiconductor microcavity [52, l23| . 

The main point of the present paper is to give a com- 
plete account of the physics of the Goldstone mode of 
a planar optical parametric oscillator above threshold, 
and to propose a way of probing it with an additional 
laser beam. In Sec^we present the physical system and 
the theoretical framework used for its description. In 
Sec IHII we study the stationary state of parametric os- 
cillation above threshold and we introduce the concept 
of spontaneous breaking of the [7(1) signal/idler phase 
symmetry. The spectrum of the elementary excitations 
around the stationary state is calculated and discussed 
in the following Sec lIVI and then used to evaluate the 
linear response of the system to a weak probe beam at 
angles close to the signal emission: the Goldstone mode 
is shown to give a strong and narrow peak at low frequen- 
cies, whose linewidth tends to zero as the direction of the 
probe beam is brought closer to the signal emission one. 
The consequences of the presence of an additional signal 
laser field at exactly the same wavevector as the signal 
emission are investigated in Sec. [V] From the full equa- 
tions of motion, one finds that the phase of the signal 
emission results in this case pinned to the one of the sig- 
nal beam so that the U(l) symmetry is explicitely broken. 



No Goldstone mode is therefore present any longer and 
a gap opens up in the dispersion of the elementary exci- 
tations. As a consequence, a dramatic broadening of the 
peak corresponding to the Goldstone mode is observed in 
the response spectrum to the probe. This phenomenol- 
ogy is the non-equilibrium analog to what happens in a 
ferromagnet when an external magnetic field is applied 
to the system to break the rotational symmetry. In this 
case, the orientation of the magnetization is fixed by the 
applied field and a gap opens up in the magnon disper- 
sion p|. 

The two last section are devoted to the analysis of is- 
sues of experimental relevance. The consequences of a 
possible frequency mismatch of the applied signal beam 
from the natural frequency of the signal emission are in- 
vestigated in Sec. I VII while the effect of the spatial in- 
homogeneities due to the finite pump spot are addressed 
in Sec lVIII Both these issues are shown not to affect the 
properties of the Goldstone mode discussed in the previ- 
ous sections. Conclusions are drawn in Sec lVlllI Three 
appendices are devoted to a brief summary of the an- 
alytical equations definining the stationary state of the 
homogeneous system, to the analytical properties of the 
Goldstone mode frequency around k = 0, and to the com- 
parison with the case of OPOs based on a second- 
order nonlinearity. 



II. PHYSICAL SYSTEM AND THEORETICAL 
MODEL 

The physical system we are considering consists of a 
planar semiconductor microcavity containing a quantum 
well with an excitonic transition strongly coupled to the 
cavity mode. A sketch is shown in Fig^ The elemen- 
tary excitations of this system are cavity-polaritons, i.e. 
coherent superpositions of cavity photons and excitons, 
which satisfy the Bose statistics. The photonic compo- 
nent is essential to couple to the external incident radia- 
tion, while the excitonic component provides the exciton- 
exciton collisional interactions which are responsible for 
the parametric process. 

Given the translational invariance of the system along 
the cavity plane, the in-plane wave vector k is a good 
quantum number and the polaritonic dispersion can be 
studied as a function of k. As shown in Fig^3, two 
polaritonic branches exist in the polariton spectrum, split 
by twice the Rabi frequency £Ir of the exciton-photon 
coupling 0,^3- A simple, but so far accurate model 
for interactions is based on a repulsive exciton-exciton 
two-body contact interaction of strength g. 

In the following we will focus our attention on the lower 
polariton branch only, which is more protected from loss 
and decoherence processes. Its dispersion will be denoted 
as e(fc). In order to justify the neglection of the upper 
polariton, one has to check that the Rabi splitting Q.r is 
much larger than both the detuning of the incident laser 
(of frequency lo p and wave vector and the nonlinear 
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FIG. 1: Upper panel: sketch of the microcavity system and 
the parametric process under consideration. Lower panel: 
Lower polariton (LP) and upper polariton (UP) dispersion. 
The cavity photon dispersion is uic(k) = ll>% \/l + fe 2 /fe| with 
hJc = 1.4 eV and k z — 20/im -1 . The exciton dispersion is 
flat and resonant with the k — cavity mode ujx = Uc- 
The exciton-photon Rabi coupling is hQ,R = 2.5 meV. Dots 
indicate the signal, pump and idler modes. 



shift e m f < griLP of the polariton mode (ulp is here 
the polariton density). These condition are actually well 
satisfied in current experiments. 

At the mean-field level, the dynamics of the polaritonic 
field can be described by a polaritonic Gross-Pitaevskii 
equation 0,0' whose A:-space form reads: 



^V'Lp(fc) = |e(fc) -i- 



+ X! 9k, qi , q2 i>* L p(qi +92 - k)ip L p(qi)ipLp(q2). (I) 

F p (k) is here the amplitude of the incident pump laser 
field, which is assumed to be continuous wave and 
monochromatic at uj p . Unless otherwise specified (as e.g. 
in Sec lVIIfl . a plane- wave at k p is taken for its spatial 
profile F p (k) = F p 5k,k ■ The damping rate 7 is the 
sum of the contribution of radiative and non-radiative 
losses. It has been taken for simplicity as momentum- 
independent. The momentum-dependent nonlinear inter- 
action strength for polaritons 3fc, gi . 92 is defined in terms 
of the Hopficld coefficient Ux quantifying the excitonic 
content of the polariton [24| as: 

9k, quQ2 = 9 U* x {k) U* x { qi +q 2 -k) U x (qi) U x {q 2 ), (2) 

Note that the same equation can be used to describe 
the photon dynamics in different systems, e.g. metal- 
lic mirror planar cavities containing a Kerr nonlinear 



medium |2jj. In this case, no excitonic field is present 
and the polaritonic field ^lp reduces to the bare e.m. 
field E. The nonlinear coupling constant g is then pro- 
vided by the Kerr optical nonlinearity of the cavity ma- 
terial (of dielectric constant eii n and thickness d), and 
is proportional to its third-order nonlinear polarizability 



X 



(3). 



% = C X (3) (^0)7(*Un), 



(3) 



where C is a numerical factor of order one that takes into 
account the boundary conditions at the cavity mirrors. 

Quantitatively, one however should keep in mind that 
the nonlinear coupling constant g for semiconductor mi- 
crocavities in the strong exciton-photon coupling regime 
hg w 5 x 10 -6 eV /im 2 is orders of magnitude larger 
than the one that is obtained in standard transparent 
media for nonlinear optics. Organic materials specifi- 
cally suited for nonlinear optical applications have in fact 
X*- 3 -* up to about 10~ 9 erg _1 c 



1 25], which yields val- 
ues for the nonlinear coupling constant of the order of 
hg = 5 x I0~ 9 eV /1m 2 . Conventional inorganic materials 
have even weaker susceptibilities. 



III. PARAMETRIC OSCILLATION AND 
SPONTANEOUS SYMMETRY BREAKING 

Among the many processes described by the motion 
equation in the different regimes, the focus of the 
present paper will be concentrated on parametric pro- 
cess in which two polaritons in the pump mode (wave 
vector k p and frequency ui p ) collide and are converted 
into a pair of signal and idler polaritons, of wave vec- 
tors respectively k s and k{ = 2k p — k s and frequencies 
uj s and u>i = 2u) p — w s . As shown in Fig^ the peculiar 
polaritonic dispersion allows for the parametric process 
to occur in a triply-resonant way, i.e. with all the lu p _ s ^ 
close to resonance with the free polariton energy e(fc PjSj i). 
Such a process is described by a term in the nonlinear in- 
teraction Hamiltonian of the form: 



Hint = gk p ,k s ,k t ^Lp(fcp) y{ P {kp) ^Lp(ki) ^ L p(k s )- 



h.c. 



(4) 



where ^ lp is the polariton quantum field operator. 

It is important to remind that OPO operation can also 
make use of a different nonlinear process, in which a sin- 
gle pump photon is split into a pair of signal/idler pho- 
tons. The nonlinear susceptibility involved in this para- 
metric downconversion process is the second-order one 
■ Before the advent of semiconductor microcavities 
in the strong coupling regime, most of the existing lit- 
erature on OPOs was indeed based on such parametric 
down conversion processes 0. Although our attention 
will be in the following concentrated on the specific case 
of semiconductor microcavities, the concept of the Gold- 
stone mode extends without almost any change to the 
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FIG. 2: Pump and signal intensity as a function of the pump 
power for k s — 0.075/im -1 . Pump parameters k p — 1.57/im^ 1 
and uj p = 1.39875. Cavity and exciton parameters correspond 
to the LP dispersion shown in figure Qd. 



other case. A brief discussion of these issues is given in 
Appendix EJ 

In order to study the parametric process the following 
ansatz can be used for ipLp(x,t): 

ij) LP {x,t) = p e H*p*-up*) +Se i(ksX-o J3 t) + 

+ Je i(fc ia - Wi t) 5 ( 5 ) 

where k s is considered here as a given parameter. A de- 
tailed analysis of its behaviour as a function of the geome- 
try and system parameters is postponed to a forthcoming 
publication [28| . 

This ansatz allows for finite amplitudes P, S, I in re- 
spectively the pump, signal and idler modes, while all 
other modes are assumed to remain empty. The possible 
occupation of other modes via multiple scattering pro- 
cesses, as observed in [2ll l29l | is therefore not taken into 
account by our model. This approximation is well justi- 
fied by the fact that the (small) polariton population in 
these modes does not affect the physics under examina- 
tion but can only give quantitatively small corrections. 
The values of P, S, I in the stationary state, as well as the 
parametric oscillation frequency u s are obtained by in- 
serting this ansatz in the equation of motion . Details 
of the equations are given in Appendix 1X1 

The pump and signal intensities in the stationary state 
are plotted in Fig|5] as a function of the pump intensity 
Ip = \F p \ 2 for a pump frequency uj p < e(fc p ), in which 
case t he p ump-only dynamics corresponds to optical lim- 
iting |18| . The main feature is the onset of parametric 
oscillation for pump intensities I p within a certain range 
of values. Outside this range, the pump-only solution 
5 = 7 = remains a dynamically stable solution of the 
mean-field dynamics. 



While the amplitude P of the pump mode is completely 
determined (both in modulus and in phase) by the inci- 
dent laser amplitude, the signal and idler phases remain 
free thanks to the invariance of the equation of motion 
and of the Hamiltonian under a simultaneous phase 
rotation of the signal and idler in opposite directions: 

S^Se lAcl> I^Ie-***. (6) 

In the parametrically oscillating state, S and I have a 
non-vanishing value, so that this £7(1) phase-rotation 
symmetry is spontaneously broken. In the absence of 
external perturbations, the specific value of the phase of 
S (and consequently of I) is randomly selected at each re- 
alization of the experiment. Consequences of the under- 
lying symmetry are however present: no restoring force 
opposes a simultaneous and opposite rotation of the sig- 
nal and idler phases, which then slowly diffuse in time 
under the effect of fluctuations [3(J. In the next section, 
we shall show how the spontaneosly broken symmetry 
manifests itself as a strong response of the system when 
spatial twists of the signal/idler phases are created by an 
extra probe laser at frequency uj r = uj s + Auj and wavevec- 
tor k r = k s + Ak close to the frequency and wavevector 
of the signal emission \Ak\ -C |fc s — k p \. 



IV. RESPONSE TO A WEAK PERTURBATION 
AND THE GOLDSTONE MODE 

Provided the applied perturbation is weak, the re- 
sponse of the polariton field can be calculated by lin- 
earising the mean-field equations of motion around the 
solution ©. Modifying this solution as 

S -» S + u s e i( ~ Akx - Au ^ +v * e -i(^kx-Ao J t) ( 7 ) 
P -> P + u p e i(Akx - Auif> + v * e - i{Akx - Aw V (8) 
I -> I + Ui e i{Akx - Au}t) +v*e- i{Akx - Au>t \ (9) 

the deviations from the steady state can be grouped in a 
6-vector U — (u s ,u p ,Ui,v s ,v p ,Vi) T that obeys the equa- 
tion 

AujU — C{Ak)U = T r (10) 

with a force vector which for our specific excitation 
scheme reads T r = (F r , 0, . . . , 0) T . The observable quan- 
tity is the number of polaritons created in the k r mode, 
which corresponds to the square modulus of the first el- 
ement u s of the system response U. 

The matrix L has the typical Bogoliubov structure 

C(Ak) - ( M{Ak) Q{Ak) ) (11) 
L - { ^> ~ \ -Q*(Ak) -M*{-Ak) )> 
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FIG. 3: Imaginary (upper panels) and real (lower panels) part 
of the excitation spectrum around the stationary state for the 
pump intensity corresponding to the dashed line in fig|5| The 
left panels refer to the case in which the signal/idler phase 
rotation symmetry is spontaneously broken and a Goldstone 
mode G is present (heavy line). The right panels refer to 
the case where a small signal laser beam {Is/I P = 0.0023) is 
present to explicitely break this symmetry and destroy the 
Goldstone mode. Cavity and pump parameters as in figure |2] 



where the 3x3 matrices M(fc) and Q(k) are defined as: 
M mn {k) = [e(k m + fc) - u m - ry/2] 5 m>n 

3 

+ 2 9k m +k,k n + k,k t Sm : n+r-tA*At (12) 

rt=l 
3 

gk m +k,k r ,k B $m+n,r+t A r A s . (13) 

rt=l 

Here we have identified the three modes s,p, i with re- 
spectively the values 1,2,3 of the matrix and vector in- 
dices, so that e.g. Ai — S, A 2 = P and A3 = I. 

The real and imaginary parts of the eigenvalues of the 
Bogoliubov matrix £(Afc) give respectively the frequency 
and linewidths of the elementary excitations around the 
parametrically oscillating stationary state. These fix the 
shape of the luminescence peaks as observed e.g. in the 
experiment of |2l| , as well as the poles of the response to 
the probe beam. A typical example of this Bogoliubov 
dispersion as a function of Ak is plotted in the left panels 
of Fig. |3| the most relevant feature is the presence of a 
branch of eigenvalues LUc(Ak) tending to exactly zero for 
Afc -> 0. 

The presence of this branch is a direct consequence 
of the fact that the mean field steady state © spon- 
taneously breaks the U(l) symmetry of the signal/idler 
phase: as any global phase rotation of S, I of the form 
© maps a stationary state of Q into another station- 
ary state with different signal/idler phases, the generator 
G T = (iS,0,—iI,—iS*,0,iI*) of the signal/idler phase 



rotations is an eigenvector of the matrix £(Ak — 0) with 
zero eigenvalue uia(Ak = 0) = 0. For finite values of Afc, 
the soft Goldstone branch at cua(Ak) corresponds to a 
spatially varying twist of the signal/idler phases. 

The existence of a soft Goldstone mode is a general re- 
sult valid for both non-equilibrium and equilibrium Q 
systems: in this latter case, the Goldstone mode corre- 
sponds to the magnon mode of ferromagnets , or to the 
zero sound mode in superfluid liquid Helium and Bose- 
Einstein condensates 0, ■ 

Fundamental differences however exist, which make 
the physics of the two cases quite distinct. In Bose sys- 
tems at equilibrium, the dispersion of the Golstone mode 
goes as Re[w(fc)] ~ c s |fc| around fc = and shows a sin- 
gularity at fc = 0. This branch corresponds to weakly 
damped pjsl ] sound waves which in the long-wavelength 
limit propagate at the speed of sound c s . 

In the present non-equilibrium case, no singularity ap- 
pears in the dispersion relation wa(Ak) of the Goldstone 
mode. In particular, the real part Re[w(Afc)] cx Afc has 
a continuous and non-vanishing slope at Afc = 0. This 
fact is due to the flow of the pump polaritons which are 
injected with a finite wave vector k p and are able to drag 
the elementary excitations. On the other hand, convec- 
tive stability and analyticity arguments (see Appendix 
lB)l show that the imaginary part goes as Im[w(Afc)] w 
—a (Ak) 2 for small Afc with a positive a > 0. 

Combining these facts together, one can physically 
summarize that the Goldstone mode of a planar OPO 
consists of a spatially slowly varying twist of the signal 
and idler phases, which however does not propagate as a 
sound wave, but has a diffusive character. Once a local- 
ized perturbation is created in the system, this will sim- 
ply relax back to the equilibrium state while it is dragged 
by the pump polariton flow. This overdamped character 
is a remarkable difference with respect to the equilibrium 
case. 

The general properties of the eigenvalues of the Bogoli- 
ubov matrix (|llfl can be used to physically understand 
the response spectra shown in FigQ] the number of po- 
laritons created by the probe in the k r — k s + Ak mode is 
here plotted as a function of its frequency uj r = lu s + Alo. 
The most apparent feature is the strong and narrow peak 
that appears at low frequencies for a probe wavevector 
k r close to the signal one k s . In particular, its linewidth 
can be much smaller than the damping rate 7 of the po- 
lariton mode: as k r approaches fc s , the linewidth goes 
to and the peak height diverges, while the peak dra- 
matically broadens for increasing Afc. Comparing these 
spectra with the dispersion of the elementary excitations 
shown in Fig|31 one can immediately see that the strong 
peak corresponds to the pole at wc(Afc), and therefore 
to the excitation of the soft Goldstone mode. 

An experimental observation of this narrow peak in 
the probe response spectra with the characteristic de- 
pendence on the probe angle would provide a unambigu- 
ous signature of the Golstone mode. In our previous 
work [23, the excitation of the Goldstone mode by the 
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FIG. 4: Response spectrum to a probe beam for several values 
of the wavevector Afc = k r — k s — 0.1, 0.15, 0.2, 0.3/^m -1 
(black, blue, red, green). Cavity and pump parameters are 
the same as in figure [5] 



quantum fluctuations showed up as a fc-space broaden- 
ing of the signal emission in the absence of any probe. 
Although this broadening of the luminescence pattern 
is of conceptual interest from the point of view of low- 
dimensional non-equilibrium physics, the pump-probe 
experiment discussed here appears as more favourable 
in view of the experimental observation of the Goldstone 
mode. The response to the probe can be measured by 
comparing pairs of spectra taken respectively in the pres- 
ence and in the absence of the probe, so to subtract out 
the background of signal light scattered by the defects of 
the sample. 

V. DESTROYING THE GOLDSTONE MODE 

It is a well-known fact that the rotational symmetry 
can be explicitely broken in a ferromagnet by adding an 
external magnetic field that fixes a preferential orienta- 
tion for the magnetization and therefore opens a gap in 
the magnon dispersion |f| . While a pinning of the Bose 
field phase is hardly obtained in Helium or atomic Bose- 
Einstein condensates because of particle number conser- 
vation [l| , it can be easily done in the optical experiment 
proposed here by applying an extra laser beam (called 
hereafter signal laser) at exactly the signal frequency u s 
and wave vector k s 39] . In this way, stimulated processes 
push the parametric emission to preferentially occur with 
a phase pinned to the incident signal laser one. This pin- 
ning of the signal field phase is the parametric and multi- 
mode analog of what happens in a single-mode laser when 
a coherent signal is injected in the cavity [32LI33T ]. In the 
pattern formation literature, the idea of forcing a spe- 
cific pattern shape by means of an additional beam often 
goes under the name of ■pattern synchronization. First 
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FIG. 5: Signal mode occupation (left panel) and amplitude 
of the gap \m[uG(Ak — 0)] (right panel) as a function of the 
signal laser power I s . The coloured points mark the I 3 values 
of the spectra in Fig[§] 




FIG. 6: Response spectrum of the polariton gas to an extra 
laser beam at Afc = 0.1 fim (left panel) and Afc = 0.2 /im 
(right panel) for several values of the signal laser intensity I a 
marked with the symbols in FiglBI 



studies of such issues in an optical context have recently 
appeared using a completely different setup |34j . 

As the Goldstone mode has a pole at exactly the signal 
laser frequency and wavevector, the change of S due to 
the application of the signal beam cannot be calculated 
by means of the linearized theory, but has to be calcu- 
lated including he signal laser beam as an additional term 
in the full equation of motion , as shown in Appendix 

El 

The left panel of Fig. [SJshows the signal amplitude as a 
function of the signal laser amplitude: for small values of 
I s , the signal intensity \S\ 2 shows a square root behaviour 
as a function of I s . As S is an analytic function of F s , 
one can expand S ~ So + cF s , the square modulus of 
which immediately gives the square-root behaviour IS^ 2 — 
\S \ 2 ^c\S \\F s \k^T s . 

The effect of the signal laser on the Bogoliubov dis- 
persion is shown in the right panels of of Fig. the 
imaginary part Imfwc^Afc)] shows now a finite gap at 
Afc = 0, so that no Goldstone mode exist any longer. 
The damping rate of Bogoliubov excitations is finite for 
any value of Afc. On the other hand the real part is much 
less affected by the signal laser beam and, in particular, 
remains at Afc = 0. 

This is another important difference as compared to 
the thermodynamical equilibrium case, where the pres- 
ence of the external field explicitely breaking the symme- 
try opens a gap in the real part of the dispersion law [j| . 
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In the present non-equilibrium case the gap appears in- 
stead in the imaginary part. The dependence of the gap 
amplitude Im[wG(Afc = 0)] on the signal laser intensity 
I s = \F S \ 2 is plotted in the right panel of Fig. [3J Again, 
the behaviour for low values of I s is a square root one. 
This can be easily explained by the fact that wc(Afc = 0) 
is an analytic function of S, P and I and these have the 
square root dependence on / s shown in the left panel of 
Fig0 

The presence of the gap in the Bogoliubov spectrum is 
expected to have dramatic consequences on the response 
of the system to the probe laser of wave vector k r and fre- 
quency uj r . This can be quantitatively studied by means 
of the same linear response equation (fTTTf) once the correc- 
tion to the stationary-state amplitudes S, P and / due 
to the presence of the signal laser have been duly taken 
into account. Note how this merely quantitative change 
in the entries of the C matrix is enough to destroy the 
weakly-damped Goldstone branch. 

Examples of probe response spectra as a function of 
the probe frequency uj r are plotted in Fig. 0for several 
values of the signal laser power / s and two different values 
of the probe wave vector k r . It is easy to see that the 
main consequence of the presence of the signal laser field 
is to broaden the peak and suppress the strongly peaked 
response to small Afc perturbations. This is in exact 
correspondance with the Bogoliubov spectra plotted in 
the right panels of FigEl 

This characteristic dependence of the probe response 
on the signal beam intensity constitutes a direct and 
experimentally accessible signature of the existence of 
a Goldstone mode which gets destroyed when the sig- 
nal/idler phases are pinned by the signal beam. 



VI. SIGNAL FREQUENCY MISMATCH 

The discussion of the previous section has assumed the 
signal laser field to be exactly on resonance with the natu- 
ral parametric oscillation frequency uj s of the cavity when 
this is illuminated by the pump laser alone. Here we shall 
study the case when the parametric emission is forced by 
the signal beam to take place at a slightly different fre- 
quency uj s + Slo s . 

Two intensity regimes have now to be distinguished for 
the signal laser. For a small intensity I s , the signal beam 
can be considered as a perturbation of the parametric 
oscillation state at the natural frequency oj s and can be 
introduced within the linear response theory of section 
IIVI as a probe beam. Thanks to the detuning 5cu s , the 
signal beam is in fact not on resonance with the Gold- 
stone mode, so that its response does no longer diverge 
as it instead happened in seclvl 

On the other hand, the linearized theory breaks down 
for larger signal laser intensity, when the parametric os- 
cillation no longer takes place at its natural frequency lo s 
but rather appears at the forcing signal beam frequency 
lu s + 5oj s . A stationary state for the mean-field equation 
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FIG. 7: The same as figure |SJ but for a signal laser detuning 
of Sujs ~ 0.01 meV. In the hatched region, the three-mode 
solution is dynamically unstable (see text). 




FIG. 8: The same as Fig. |SJ but for a signal laser detuning 
of Scus = O.OlmeF and for the values of 7 S indicated with the 
symbols in Fig'0 



JU of the form: 
i/> LP (x, t)=P e <(M-<"**) + 5- e i{k,x-{u>,+Su,)t) + 

_|_ j e i(kiX-(ui-&u s )t) ^ ^4) 

has thus to be considered. The transition from one case 
to the other can be described by means of a Hopf bi- 
furcation scenario in a theory explicitely including the 
possibility of a time-dependence for S, P and I. For the 
sake of simplicity, we will restrict ourselves in the fol- 
lowing to the parameter range well above the bifurcation 
point, where the ansatz 1|14[) is accurate. The Hopf bifur- 
cation point is signalled by the solution (|14fl becoming 
unstable, i.e. the imaginary part of an eigenvalue of the 
linearized theory becoming positive Im[w(Afc = 0)] > 0. 

The signal mode occupation n s and the imaginary part 
Im[o;G(Afc = 0)] of the Goldstone mode at Afc = are 
plotted in respectively the left and right panels of Fig. 
as a function of the signal beam intensity I s . The 
hatched region is the low I s region where the solution 
(|14fl is dynamically unstable. Well above the instability 
region, the behaviour of Im[ujQ(Ak = 0)] as a function 
of the signal laser intensity recovers the behaviour of the 
zero-detuning 5u> s — case. 

Examples of response spectra to the probe beam are 
shown in Fig|S] For low values of the signal laser in- 
tensity I s (but still above the instability threshold) the 
response is strongly peaked at the weakly damped Gold- 
stone mode, while it broadens for larger values of I s : this 
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FIG. 9: Numerical results for the spatial profile of the signal 
density (left panel) and the signal wave vector spectrum (right 
panel) for a pump waist a v — 300 /im, a carrier pump wave 
vector k p — 1.57 /im. The vertical dashed line in the right 
panel shows the wavevector of the signal emission at the center 
of the spot. Other parameters as in Fig|5] 



phenomenology is in close qualitative agreement with the 
one shown in Fig. for the case of a perfectly tuned sig- 
nal laser. Provided we are sufficiently far from the Hopf 
bifurcation, a mismatch in the signal frequency therefore 
does not substantially affect the signal/idler phase pin- 
ning effect which is responsible for the disappearance of 
the strong response associated to the Goldstone mode. 
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(co r -co s )/y 

FIG. 10: Solid line: response of the polariton pattern of Fig. 
[§|as a function of frequency, when it is probed by a laser beam 
with wave vector k r — 0.175/im _1 (full black line). Dashed 
line: prediction for a homogeneous system with the values at 
the center of the spot. Red dotted line shows the response 
when a signal laser with intensity I B = 10 -3 I p and waist 
o~s ^> o v is applied to the inhomogeneous system. 



VII. FINITE SPOT 



In all the previous discussion, a spatially homogeneous 
system was considered, with a spatially homogeneous 
pump beam; in that case, the wave vector was a good 
quantum number. Real experiments are however per- 
formed using finite-size pump laser spots, generally of 
Gaussian shape F p (k) ~ exp[— a p (k — k p ) 2 /2]. This im- 
plies that a range of wave vectors are excited, yet the 
pump, signal and idler beams remain well distinct in 
k space provided the (real-space) spot size a p is large 
enough \k p — k g \ a p 3> 1. Under such an hypothesis, the 
signal and idler phase rotation © is still a symmetry 
element of the problem j4jj • 

The spatial and wave vector distributions of the para- 
metric signal emission are shown in Fig. [!|]for the case of 
a wide pump spot in the absence of any incident signal 
laser. The plotted curves are the result of a numerical 
integration of the full mean-field evolution equation i[TJl 
until the steady-state is reached: the main feature to 
note is the inhomogeneous wave vector broadening that 
can be observed in the right panel, a broadening which 
always remains much smaller than the wavevector space 
distance \k p — k s \ ~ 1.5 /im" 1 . A detailed explanation of 
the physical origin of this broadening is posponed to the 
forthcoming publication psj . 

The response of this steady-state to an additional 
probe laser beam with wave vector k r = 0.175 and a waist 
oy = cr p is shown in Fig. ^] As in the case of an homo- 
geneous system, the very peaked response as a function 
of frequency is a consequence of the presence of the Gold- 
stone mode. The enhanced broadening as compared to 



the homogeneous case is mostly due the inhomogeneous 
broadening of the signal emission in wavevector space 
shown in figEl The top of the peak is in good agreement 
with the prediction for a homogeneous system with the 
same signal wave vector as the one actually selected at 
the center of the finite spot (dashed line). Agreement 
is observed also in the low-frequency tail of the peak, 
but differences are apparent in the high-frequency tail: 
the tail of the finite spot curve extends in fact for much 
longer than in the homogeneous case. This is a direct 
consequence of the asymmetric inhomogeneous broaden- 
ing of the signal wave vectors shown in the right panel 
of Fig. [51 which indeed mostly extends towards lower k s 
values. Relative to these wave vectors, the perturbation 
has a larger Ak = k r — k s , which means that the response 
peak lies at higher frequencies and has a larger width (as 
shown in the lower panels of Fig[2J). 

The red dotted line in Fig. shows the response 
when a signal laser beam is applied, which explicitely 
breaks the U(l) symmetry. Analogously to the case of a 
homogeneous system, the fact that the U(l) symmetry 
is now explicitely broken causes a strong broadening of 
the peak and a corresponding suppression of the peak 
response. This completes the proof that the concepts 
and the analytical results developed for the homogeneous 
case can still be safely applied when the pump spot has 
a finite size. 
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VIII. CONCLUSIONS 

In the present paper, we have studied the properties of 
the Goldstone mode that appears in the elementary ex- 
citation spectrum of planar parametric oscillators above 
threshold as a consequence of the spontaneous breaking 
of the U(l) signal/idler phase symmetry. Although quan- 
titative predictions have been shown only for the case of 
planar semiconductor microcavities with a quantum well 
excitonic transition strongly coupled to the cavity mode, 
the physics is the same in any kind of parametric oscilla- 
tor with planar geometry. 

Analogies and differences with thermodymical equilib- 
rium systems showing the same spontaneous symmetry 
breaking phenomenon (e.g. superfluid Helium 4 or Bose- 
Einstein condensates) have been drawn. In particular, 
several features have been pointed out, which are typi- 
cal of the stationary state of a non-equilibrium system 
where a dynamical equilibrium occurs between driving 
and losses. Although its dispersion tends as usual to 
in the long wavelength limit, the Goldstone mode of our 
non-equilibrium system is not a propagating mode like 
zero-sound in Helium or Bose-Einstein condensates, but 
rather an overdamped mode. A simple way of probing 
its dispersion at small wavevectors is to measure the re- 
sponse of the system to an additional probe laser shone 
onto the cavity at an angle close to the signal emission. 
As the probe beam approaches the signal, the damping 
rate of the Goldstone mode tends to vanish, which corre- 
sponds to a probe response spectrum extremely peaked 
at the Goldstone mode frequency. 

A signal laser beam injected into the cavity and stim- 
ulating the parametric emission is shown to provide an 
external field which explicitely breaks the U(l) symme- 
try. This in analogy with what happens in a ferromagnet, 
where an external magnetic field provides a preferential 
orientation to the magnetization. As usual, the explicit 
breaking of the symmetry is responsible for the opening 
of a gap in the Goldstone mode dispersion. Differently 
from the equilibrium case, the gap is here in the imagi- 
nary part of the dispersion. As a consequence, the peak 



in the response spectrum gets broader and eventually al- 
most disappears when the signal laser intensity grows 
higher. 

In the last sections of the paper, a few issues of experi- 
mental relevance in view of the actual observation of the 
Goldstone mode have been addressed, such as the effect 
of a slight detuning of the signal laser from the natu- 
ral parametric oscillation frequency, and the finite size of 
the excitation spot. All of them are shown not to give 
any dramatic effect on the observation of the Goldstone 
mode. 

An experimental study along these lines appears there- 
fore feasible with present technology samples and will 
constitute an important step in the understanding of 
the properties of parametric oscillation in spatially dis- 
tributed systems, in particular of its analogies and differ- 
ences with the Bose-Einstein condensation phase transi- 
tions in systems at thermodynamical equilibrium. 
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APPENDIX A: EQUATIONS DEFINING THE 
STATIONARY-STATE 

Substituting the ansatz (J5J in the equations of motion 
JTJ and imposing the stationary state leads to the follow- 
ing set of complex equations: 
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P + (\P\ 2 + 2\S\ 2 + 2|/| 2 ) P + 2P*SI + F p 
S + (2\P\ 2 + \S\ 2 + 2|/| 2 ) S + P 2 I* 

1+ (2\P\ 2 + 2\S\ 2 + |/| 2 ) I + P 2 S*. 



(Al) 
(A2) 
(A3) 



r 



The following shorthand notations have been introduced F p — yfg Uc(k p ) /Ux (k p ) F p have been defined. The com- 
ep,s,i — £{k P ,s,i) and X PtS ,i — Ux(k p . s ,i)- Scaled quan- pi ex variables S, P, I, u s can be obtained from the set of 
tities S = ^/gX s S, P = ^/gX p P, I = ^fgXil and three complex equations (|A1IA3|) together with the con- 
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dition that the frequency lu s has to be real: only 7 real 
equations are then available to determine 8 real variables, 
which means that the solutions are grouped in ID man- 
ifolds. As the set of equations f l|All) - (|A3|) ') is symmetric 
under the phase rotation JfjJ), the solution manifold is 
generated by the action of © on a given solution. 

The presence of the signal beam at frequency uj s + 8w s 
is taken into account in the set of equations (' ljAl|) - (|A3|) ') 
by simply adding the term F s = ^fgUc{k s )/Ux{k s ) F a to 
the right-hand side of (|A2() and considering lo s — » u> s ±5u) s 
as a fixed quantity. In this case, the number of unknowns 
is reduced to three complex variables, which are then 
completely determined by the three complex equations. 



APPENDIX B: ANALYTICAL REMARKS ON 
THE GOLDSTONE MODE 

In this appendix we show that (convective) stability |6J 
implies that only even powers are possible in the expan- 
sion of Im[w(fc)] around k = 0. The eigenvalues of a lin- 
ear operator C{k) have a Taylor expansion around every 
value of k, except for some spec ial p oints where fractional 
powers 1/p can occur (see Ref.[36j], p. 65). If the Taylor 
expansion exists, stability requires Im[w(fc)] < for all 
k, which implies that only even powers of k are allowed 
in the Taylor expansion of Im[w(fc)] around k = 0. 

If accidentally k — is a special point of order p, there 
are p eigenvalues which have a Puiseux series instead of 
a Taylor series 



where a single pump photon is converted into a pair of 
signal/idler photons. An ansatz analogous to © can 
then be used for the stationary state of the three pump, 
signal and idler fields at mean-field regime: 



E p (x) 
E s (x) 
Et(x) 



p ^i(kpX-ujpt) 
j i{kiX—uJit) 



(C2) 
(C3) 
(C4) 



and the amplitudes P, S, and /, as well as the frequency 
uj s are obtained by inserting this ansatz into the field 
equations of motion for the Hamiltonian ijClfl . 

It is appearent that the transformation © is still a 
symmetry of the Hamiltonian (|C1|) . which is sponta- 
neously broken above threshold by the ansatz (|C2IC4(1 . 
A Goldstone mode is therefore present, with the same 
properties as discussed in the body of the present paper. 

The only different lies in the dependence of the sta- 
tionary state on the pump intensity |37|: as the second- 
order nonlinearity does not produce any direct shift of the 
mode frequency due to the mode population, the curves 
of FigI3 are dramatically modified. In particular, there 
is no upper threshold for the parametric process, which 
at high intensities is replaced by chaotic behaviour [37|. 

In the degenerate case where signal and idler belong 
to the same photonic L branch, the parametric terms of 
the parametric Hamiltonian read instead: 



hg 



M{k v ) E L {h) E L (k s ) +h.c. 



(C5) 



X h (k) = A + a 1 e 27Ilh / p k 1 / p + a 2 e 4mh/p k 2 / p + • • • , (Bl) so that the ansatz 



for h = 0, 1, . . . ,p — 1. If p = 2, stability requires that 
Im(a m ) = for any odd m, so that no half- integer powers 
are allowed in the expansion of Im[A(w)]. The case with 
p > 2 is ruled out just because it would forcedly give an 
eigenvalue with positive imaginary part. 

In summary, only even powers are allowed in the ex- 
pansion of Im[a;(fc)], so that the expansion of Im[w(/;)] is 
to leading order quadratic (unless this term accidentally 
vanishes). 



APPENDIX C: COMPARISONS WITH x (2) OPO'S 

In this Appendix we make the link of our theory to the 
different case of an optical parametric oscillator based 
on a nonlinear medium showing a second-order optical 
nonlinearity instead of a third order one |S( • 

In the non-degenerate case where pump, signal and 
idler belong to different photonic branches of the planar 
cavity, the terms of the Hamiltonian responsible for the 
parametric downconversion process have the form: 



H,, 



tig 



mk p )Ei(ki)E s (k s )+b.c. 



(CI) 



%{x) -» Pe** -1 *') (C6) 
ip L (x) -» s e i(*.3-<".t) + j e «(*i*-"it) (C7) 

has to be used to describe the parametric oscillation 
state. Also in this case, the equations of motion are in- 
variant with respect to the transformation © and all the 
physics of the Goldstone mode remains unchanged. Note 
how the transformation (|HJ corresponds in this case to a 
spatial translation of the field (|C7jl in the L mode. 

The only exception occurs when the system oscillates 
in a completely degenerate regime k s = ki = k p /2 where 
the signal and the idler coincide: in this case, the U(l) 
symmetry is reduced to a discrete ± symmetry and no 
Goldstone mode is any longer present. A linear analy- 
sis for this completely degenerate case has been worked 
out in 37], where it was shown explicitely that no zero 
eigenvalue is present for Afc = 0. 

Note that the condition characterizing a complete de- 
generacy regime is weakened for a finite spot of size cr p , 
where the pump, signal and idler spots have a /c-space 
broadening of the order of l/a p . In this case, in fact, the 
Goldstone mode can be shown to disappear as soon as 

I k s k>i | — 1 . 
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